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ABSTRACT. A novel representation of a quasi-periodic modified von Mangoldt 
function L(n) on prime numbers and its decomposition into Fourier series has been 
investigated. We focus on some particular quantities characterizing the modified 
von Mangoldt function. The results indicate that prime number progression can be 
decomposed into periodic sequences. The main approach is to decompose it into 
sin or cosine function. Basically, it is applied to extract hidden periodicities in 
seemingly quasi periodic prime function. Numerical evidences were provided to 
confirm the periodic distribution of primes. 
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I. The von Mangoldt function 


The von Mangoldt arithmetic function is defined by 


h(ji) — 0°g V if tl — p k for p a prime and 
l 0 otherwise. 


iwh/i n n'v 
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This function carries the information about the properties of the primes and 
approximate their weighted characteristic function. The discrete Fourier transform 
of the von Mangoldt function is defined by the Riemann-Weil explicit formula 


lr m = h (j) + h(-j)- 

—^(O) log n ++y Mr+(++)*- - 2 XC.+aflogn) 2 

and gives a spectrum with spikes at ordinates equal to the imaginary part of the 
Rienrann zeta function zeros (Hejhal 1976). 
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II. The modified von Mangoldt function 

Starting with the function above, let us consider a special, simple case of von 
Mangoldt function at for k — 1, 

L(ri) = [ l n ( e ) if n — p for p aprime 
tO otherwise. 

The function L (n) is a representative way to introduce the set of primes with 
weight of unity attached to the location of a prime. From the summation of that 
function we can see that primes contribute equally and that representation exist. 
The sum of the modified von Mangoldt function L(n ) is same as the Riemann 
hypothesis (Devenport 1980, Vardi 1991) with a slight error: 

*P(n) = / 2 P ”dt = n + 0(Vn (In x) 2 ) = S n ^oob(n) = p n 4 


but the L(n) function predict the number of primes with no error, e.g. the primes 
are not randomly distributed (large numbers’ law). Therefore, this function is 
oscillatory but not diverging. The modified von Mangoldt function L(n) is related 
to Dirichlet seriesor the Riemann zeta function by 


(C a ) = Zn =1 
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L(n) 

n a 


91(a) > 1. 
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The discrete Fourier transform of the modified von Mangoldt function L(n ) gives 
a spectrum with periodic, real parts as spikes at the frequency axis ordinates. 
Hence the modified von Mangoldt function L(n) and the prime series can be 
approximated by periodic sin or cos functions. Let us consider that natural numbers 
form a discrete function g (n) —> g(n k ) and g k = g(n k ) by uniform sampling at 
every sampling point n k = kA, where A is the sampling period and k = 2, 3,.., L — 
1. If the subset of discrete prime numbers is generated by the uniform sampling 
of a wider discrete set, it is possible to obtain the Discrete Fourier Transform 
(DFT) of L(n) . The sequence of L(n) prime function is transformed into another 
sequence of N complex numbers according to the following equation: 

F{L(n)} = X(v ) 6 


where L(n) is the modified von Mangoldt function. The the operator F is defined 
as: 


*(v) = 7 

where v — lAf, l = 0,1,2,... ,N — 1, and Af—l/L, L=length of the modified von 
Mangoldt function. 


2 



The amplitude spectrum of equation 7, describes completely the discrete-time 
Fourier transform of the L(n) prime function of quasi periodic sequences, which 
comprises only discrete frequency component. The important properties of the 
amplitude spectrum are as follows: 


1) the ratio of consecutive frequencies 1 ) and maximum frequency ( f max ) 

ft +1 ft' Kl/tO fmax ■* , r 1 , ,i 

converges lim -—-—-— 1, where f s — - represents the 


t—>f s /2 ft ft ft 

positive half-length of amplitude spectrum. This is similar to the prime 


number theorem: 


7r(x)lnx 

lim — = 1. 
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2) the reciprocate of the frequencies ( f t ) converges lirn — = — 9 

3) the consecutive frequencies on the DFT amplitude spectrum describe a 

parabolic or Fermat’s spiral (Fig.l.), starting from 0 with increasing modulus 
(polar form r = af t ), where a — 1, independently from the length of the 
prime function and its parametric equations can be written as: x(/ t ) = 
ft cos (f t ), y(f t ) = f t sin(f t ) 10 



Fig. 1. Parabolic or Fermat spiral using the consecutive frequencies. The arc length 
of the spiral is a function of the sampling interval. 

4) the DFT spectrum of the primes exhibits highly ordered and symmetrical 
frequency distribution (Fig.3.) 

5) the DFT of the modified von Mangoldt function will be periodic (Fig.2.) 
independently the length of integer sequences used (v > N,v — 0,1, ...,N — 1), 

X(v ) = X(v + IV) = X(v + 2N),... , = X(v + zN ) 

X(zN + v) = £n=o L(n) e -* (27mvXziV+v) = Yn=o Lfrfie-^^e- 1271271 11 

z and n are integers and e~ l2nzn = 1. 

X(v + zN) = Y.n= o7(n)e“ i(27rnv) = X(v) 12 

III. Spectral evidence 

For completeness we give the compelling spectral evidence for the truth that 
Fourier decomposition of the modified von Mangoldt function is periodic. Fourier 
transformation has been used for studying the behaviour of prime numbers for a 
long time and it was employed here to reveal the periodic distribution of the 
primes. Figure 2 shows that the corresponding amplitude spectrum of the modified 
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von Mangoldt function exhibits a well-defined periodic signature. This suggests that 
the series of prime numbers is not an arbitrary function. Creating the superposition 
from the Fermat spiral frequencies by adding their amplitude and phase together, 
the original prime function can be restored. 



Fig. 2. Representation of the modified von Mangoldt function up to 100. 
Spikes appear at the location of prime numbers. 



F requency 

Fig. 3. The Discrete Fourier amplitude spectrum of the modified von 
Mangoldt function a) DFT of primes between 2-10.000 (continuous black 
line) and b) DFT of primes between 1.0200*10''7-1.0201*10''7 (dotted grey 
line) 
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Fig. 4. Reconstruction of the modified von Mangoldt function from the 
Fourier decomposed periodic waves between 2-10.000 (magnified from 2-30). 
The circles at the top of the spikes indicates the location of prime numbers. 


Creating the superposition from the Fermat spiral frequencies by adding their 
amplitude and phase together, the original prime function can be arising. 

L (n) = Xsin(A + 13 

From the wave sequence above, it can be seen that all primes have a last digit of 
1, 3, 5 and 7 or 9 [4], so the music of the primes billow over the set of natural 
numbers and able to show some arithmetical progression, pattern of prime numbers 
[5]. Figure 5 shows that new prime numbers are decreasing to infinity and that 
how many new primes we have for example in a 1000 number around 1000, 10"3, 
10'4 etc., so less and less amount of new frequencies requires to describe their 
pattern. 


5 




Fig. 5. Number of new primes in a 1000 number interval. 


Conclusion 

We have presented a novel method for the detection of prime repetitions on the 
set of natural numbers. The method is based on a DFT of the modified von 
Mangoldt function. The relevance of the method for determination of the 
conformational sequences is described. 
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